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Abstract 

In this work we study the polytope associated with a 0/1 integer programming 
formulation for the Equitable Coloring Problem. We find several families of valid 
inequalities and derive sufficient conditions in order to be facet-defining inequali- 
ties. We also present computational evidence of the effectiveness of including these 
inequalities as cuts in a Branch & Cut algorithm. 
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1 Introduction and preliminary results 

The Equitable Coloring Problem (ECP), originally presented in [2], is a vari- 
ation of the widely studied Graph Coloring Problem (GCP) with additional 
constraints imposing that any pair of color classes has to differ in size by at 
most one. Further references and applications can be seen in [1]. 
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A k-coloring of a graph G = (V, E) is a partition of V in k stable sets, 
Cj, with 1 < j < k. The stable set Cj is the class of color j. An equitable k- 
coloring (or just fc-eqcol) of G is a /c-coloring satisfying the equity constraints, 
i.e. |_ n /^J < | j < \n/k] for each 1 < j < k, where n = \V\. 

Unlike GCP, a graph admiting a fc-eqcol may not admit a (fc + l)-eqcol. 
This leads us to define (G) as the set of k < n such that G does not admit 
any /c-eqcol. For instance, ^(i"^) = {1,3}. 

The equitable chromatic number of G, Xeq{G), is the minimum k for which 
G has a /c-eqcol. Computing x eg (G) for arbitrary graphs is an iVP-hard prob- 
lem [1]. 

Although many integer programming formulations are known for GCP, as 
far as we know, just two of these models were adapted for ECP. One case is 
the model in [3], adapted in [5]. Preliminary results concerning a Branch & 
Cut algorithm based on one of the models in [4] were presented in [6]. The 
algorithm turns out to be competitive compared to the one presented in [5]. 
This encouraged us to delve into a polyhedral study with the aim of finding 
strong inequalities that allow us to improve the performance of our algorithm. 

2 The polytope ECV 

From now on, we assume that G is a graph with n vertices such that n > 5 
and 2 < Xeq{G) < n — 2. Other cases are trivial. 

In [4], colorings of G are identified with binary vectors (x,w) G {0, l} n2+n 
where x G {0, l} n and w G {0, l} n , satisfying the following constraints: 

n 

x V j = 1 V v G V (assign a unique color to each vertex) 

i=i 

x u j + x v j < Wj V uv G E, j = 1, . . . , n (adjacent vertices do not share the same color) 
w j+i < w j V j = 1, . . . , ft — 1. (eliminate some symmetric colorings) 

where x V j = 1 if color j is assigned to vertex v and Wj = 1 if color j is used, i.e. 
Cj 7^ 0. The coloring polytope CP is defined as the convex hull of colorings of 
G. In this work, equitable colorings are identified with binary vectors defining 
colorings which also satisfy 

x V j <Wj, Vd isolated, j = 1, . . . , n, (1) 



k=j 



v£V k=j 



(wk-Wk+i), Vj = l,...,n- 1, (2) 



where w n+ i is a dummy variable set to 0, constraints (1) ensure that isolated 
vertices use enabled colors and (2) are the equity constraints. The Equitable 
Coloring Polytope ECV is the convex hull of the equitable colorings of G. 



Next we state the main results related to the polyhedral structure of SCV. 

Proposition 2.1 The dimension of SCV is n 2 — (\jtf(G) \ + 2). 

In [4], clique inequalities and block inequalities are proven to be facet- 
defining inequalities of CV. In our case, we have: 

Proposition 2.2 (i) Let j < n — 1 and Q be maximal clique of G such that 
\Q\ > 2. Then, the clique inequality Yliv&Q Xr "j — w j defines a facet of SCV. 
(ii) Let v G V and j < n — 2. Then, the block inequality Ylk=i 
valid for SCV and defines a facet of SCV if j — 1 ^ srf (G). 

By lifting rank inequalities and neighborhood inequalities, also studied in 
[4], we obtain new families of valid inequalities which often define facets. 

Proposition 2.3 Let j < n — 1, S C V with a(S) = 2 and Q = {q : q E 

S, S C iV[<7]}. Then, the (S, Q)-2-rank inequality defined as 

Yl Xvj+t^Xvj < 2 Wj . 

v£S\Q v£Q 

is valid for SCV. Let us assume that \Q\ > 2 and no connected component of 
the complement graph of G[S\Q] is bipartite. The inequality defines a facet of 
SCV if one of the following conditions holds: 

• for all v G V\S, Q U {v} is not a clique, 

• n is odd, j < \n/2] — 1 and for all v G V\S such that Q C N(v), there 
exists a stable set H of size 3 such that v G H and \H D S\ = 2, and the 
complement of G — H has a perfect matching, 

• n is even, j < \n/2\ — 1 and for all v G V\S such that Q C N(v), there 
exist two disjoint stable sets of size 3, H and H' , such that v G H and 
\H fl S\ = 2, and the complement of G — (H U H') has a perfect matching. 

If Q = or Q = {q}, the (S, Q)-2-rank inequality is respectively domi- 
nated by the inequalities 

^ x V j + ^ x vn -i < 2Wj + w n ^i - w n , or 

^ ^ Xyj ~\~ *2Xqj ~\~ Xqn ^ *2.Wj 
v&S\{q} 

which also usually define facets of SCV. 

Proposition 2.4 Given j < n — 1, u G V and S C N(u) with a(S) > 2, the 
(u, j, S^-subneighborhood inequality defined as 



Usxuj + ^2 Xv i + ^2 (^js - lks)x u k < Ijswj, 

vSS k=j+l 

where j k s — min{ \n/k] , a(S)}, is a valid inequality for ECV. If S = N(u) 
or a(S) < \n/j~\ — 1, the inequality defines a facet of ECV when the following 
conditions hold: 

• for all k G {["] — 1 : 2 < i < 7^-5 — 1}, there exists a k-eqcol such that 
\Cj nS\= 7fc 5 , 

• for all v G N(u)\S, there exists an equitable coloring such that \Cj fl S\ = 
a{S) and (Cj n N(u))\S = {v}. 

Finally, we obtain three new families of valid inequalities for £CV, which 
were not derived from any of the valid inequalities given in [4]. 

Proposition 2.5 Let S C {1, . . . ,n}. The S-color inequality defined as 

n 

y~l - X] b sk{wk - w k+ i), 

jeSvev k=i 

where d Sk = \S H and h Sk = d Sk [j\ + mm{d Sk ,n - k[^\}, is a 

valid inequality for ECV. In addition, if 3 < \S\ < n — 2, S contains all the 
colors greater than n — and the complement of G has a matching of 

size \ 2 "1; then the S -color inequality defines a facet of ECV. 

Proposition 2.6 Given u a non universal vertex of G and j < \n/2\ such 
that a(N(u)) > \n/j\, the (u,j)-outside-neighborhood inequality defined as 

n n 

([n/j\ - l)x uj - ^ x vj + b jk x uk< ^2 b jk(w k - w k+i ), 

v£V\N[u] k=j+l k=j+l 

where bj k = \n/j\ — \n/k\, is valid for ECV and defines a facet of ECV if the 
following conditions hold: 

• there exists v G V\iV[w] such that N(u)\N(v) ^ 0, 

• if n is odd, the complement of G — u has a perfect matching, 

• for all v G V^\A r [w], there exists a \n/2\-eqcol such that Cj = {u,v}, 

• for all k such that j < k < \n/2\ and |_fj > |_FPi J > there exists a k- 
eqcol such that \Cj fl N(u)\ = \n/k], and a k-eqcol such that u G Cj and 
\C,\N[u]\ = [n/k\-l, 

• for all k G {j, . . . ,n — 3}\^/(G), there exists a k-eqcol lying on the face 
defined by the inequality. 



Proposition 2.7 Given u G V , Q be a clique of G such that Q n N[u] = 
and j, k such that j ' < k < n — 2 and a(N(u)) > \n/k] — 1. The (u,j, k, Q)- 
clique-neighborhood inequality defined as 

n 

(\n/k]-l)x uj + x vj + ^ (\n/k]-\n/l])x u i + y^x vn _i+ ^ x vn 

veN(u)UQ l=k+l vGV veV\{u} 

k—1 n—2 n 

< ^ Kljwi - u>i +1 ) + ^ [n/fc] (w; - wi +1 ) + y^ ([n/^l+l)^-^!), 

i=j i=fc l=n— 1 

where b u \ = min{\ n/ 1], a (N(u)) + 1}, is a va/zd inequality for SCV. If there 
exists v G Q sitc/i i/iai N(u)\N(v) ^ 0, the inequality defines a facet of SCV 
when the following conditions hold: 

• for all I G {j, . . . , n — 3}\£/(G), there exists an l-eqcol lying on the face 
defined by the inequality, 

• for all v G V\(A r [w] U Q), there exist two k-eqcols lying on the face defined 
by the inequality, with v G Cj in the first one and where the second one is 
obtained from the first by only changing the color of v, i.e. v ^ Cj, 

• for all 1 < i < \n/k~\ —1, if I = max{[y] — l,n — 2}, there exist two l-eqcols 
lying on the face defined by the inequality such that u G Cj in one of them 
and u G C\ in the other. 

Although the sufficient conditions in the previous results are strong, we 
find several cases where they hold. Moreover, even when the inequalities 
do not define facets, the dimension of the faces defined by them is quite 
high. For example, if k < [n/2] — 1, it can be proved that the dimension 
of the face defined by the (u, j, k, Q)-clique-neighborhood inequality is at least 
dim(SCV) - (3n - \^{G)\ - [n/2\ - \N(u)\ - \Q\ - 5). 

3 Computational performance of valid inequalities 

In this section, we report on the computational performance of the families of 
valid inequalities studied in the previous section, embedded as cuts in a B&C 
algorithm for solving ECP. 

In order to strengthen the formulation and avoid considering classes of 
symmetric colorings, constraints x V j = 0, V 1 < w < j < n are considered 
within the initial relaxation, and x V j < Ylu=j-i x, uj-ii V2<j<t/<n are 
handled as cuts during the optimization. 

The cutting process consists in looking for violated clique and (S, Q)-2-rank 
inequalities with a greedy algorithm. During the separation of clique inequal- 
ities, it attempts to find violated (u, j, k, Q)-clique-neighborhood inequalities 



by scanning vertices u not adjacent to a given clique Q. Whenever not enough 
cuts were generated, it tries to add block, (u, j, A r (w))-subneighborhood and 
{u, j)-outside-neighborhood inequalities, handled by enumeration, and S'-color 
inequalities with a greedy algorithm. Separation routines for clique and block 
inequalities are exposed in [4]. The B&C algorithm also includes an initial 
heuristic, a primal heuristic and a custom branching rule. 

Experiments were carried out over random instances of 70 vertices with 
different density percentages and 2 hours time limit. We compare our B&C 
algorithm with (BC + ) and without (BC) our new inequalities against the 
general purpose IP-solver CPLEX 12.1 and results reported in [5]. 



% 


% solved inst. 


Nodes (average) 


Time in 


sec. (average) 


dens. 


BC+ BC CPX [5] 


BC+ BC 


CPX [5] 


BC+ BC 


CPX [5] 


10 


100 100 100 100 


3.4 4 


13.3 57 


0.3 0.3 


4 109 


30 


90 90 


2135 3949 




276 224 




50 


70 70 


7932 21595 




1354 2145 




70 


80 80 10 100 


525 2970 


214 678 


128 446 


4380 273 


90 


100 100 100 100 


5.1 14.5 


30 9.4 


2.6 2.8 


29 11 



As one may appreciate from the table, the addition of our cutting planes 
has shown to be particularly useful in substantially decreasing the number 
of Branch-and-Bound nodes and the CPU time was significantly reduced on 
medium and high density instances. 
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